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ABSTRACT 


This report summarizes the work performed duirog the third six- 
month period of the NASA Lewis sponsored research program (Contract 
Number NAG 3-67) entitled, "Prediction of Sound Radiation from Different 
Practical 3et Engine Inlets". The work done during tlie first year of the 
contract may be summarized as follows. During the initial six-month period 
the computer codes necessary for this study were developed and checked 
against exact solutions generated by the point source method using the 
NASA Lewis QCSEE inlet geometry. During the second half year, these 
computer codes were used to predict the acoustic properties of the following 
five inlet configurations: the NASA Langley Bellrnouth, the NASA Lewis 
3T15D-1 Ground Test Nacelle, and three finite hyperbolic inlets of 50, 70 
and 90 degrees. 

During the past six months, 35 computer runs were done for the 
NASA l.angley Bellrnouth. For each of these computer runs, the reflection 
coefficient at the duct exit plane was calculated as was the far field 
radiation pattern. These results are presented in both graphical and tabular 
form in this report with many of the results cross plotted so that trends in 
the results verses cut-off ratio (wave number) and tangential mode number 
may be easily Identified. 



Table of Contents 


Pace 

iiiim 1 fi Iril ■! 

Introduction I 

Numerical Considerations 3 

The Computer Runs 

Reflection Coefficients 6 

Radiation Patterns 8 

References 10 

Tables 12 

Figures I 


Introduction 


During the first year of this contract, the two developmental tasks 

associated with this project were completed. Firstly, existing computer 

codes for the calcula'.ion of the sound radiated from axisymnr/etric bodies 

were upgraded and refined so that calculations could be made having non- 

dimensional wave numbers (based on duct radius) of up to 20. This initial 

development work was accomplished using the NASA Lewis QCSEE inlet of 

Ref. 1. The results of this initial development work are presented in the 

NASA Lewis Semi-Annual Status Report (6/2/80-11/20/80) for this grant 

(Number NAG 3-67). Secondly, many small satellite programs were 

developed to do such things as calculate the cut-off wave numbers fo 

various modes, calculate the modal source inputs for the other computer 

programs, and plot the data generated in the field. These and the two main 

programs, for calculating the surface distributions and field distributions of 

the acoustic quantities of interest, were then combined into a coherent 

package of programs (i.e. the data file structure for each is the same so that 

only a single set of data files is required). This package of programs was 

then used to investigate the acoustic behavior of five different inlet 

2 3 

configurations, namelyt the NASA Langley Bellmouth ’ , the NASA Lewis 

u 5 

.1T15D 1 ground test nacelle ’ , and three finite hyperbolic inlets similar to 

those of the Y CHO^’^ study of 50, 70 and 90 degrees. Comparisons 

8 “ 1 0 

were made between the results generated by the integral technique 
used in this study and those of other experimental and theoretical studies for 
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the five inlet confij?,urations. Both reflection coefficients and radiation pat- 
terns in the field were compared for various modal inputs at various wave 
numbers. These results are presented in the NASA Lewis Semi-Annual 

Status Report (12/1/80-5/10/81) tor this grant. 

[luring the past six months, 15 computer runs were done for tlie 

NASA Langley Bellmouth. These runs form a parametric study for this 

geometry which c'onsists of five different input modes vind six cut-off ratios 

plus five special runs at other cut-off ratios done at tlie sponsors request. 

For each run, the reflection coefficient at the duct exit plane is cali'ualt 'd 

along with the c^istribution of the acoustic sound power level (SPL) in the 

field. These results are plotted and cross plotted holding either tht input 

mode or the cut-off ratio constant. 
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Numerical Con!»iderations 


Since the required computer runs for the parametric study for this 
inlet geometry cover a wide range of input modes and cut-off ratios (i.e. 
wave numbers) the same number of calculativtnal points on the body was not 
used for all of the runs. In Table I is found the range of values of tlie input 
parameters for which the efferent numbers of points vtlong the body and in 
the 0 (tangential) dl ection are used. As the wave number (cut-off ratio) 
is Increased, more points must be used on the body to accurately describe 
the distribution of the acoustic qikantltles on the surface of the body as is 
Indicated in Table I. Also, as the input mode number (M) is increased, the 
number of points used In the O integration must be increased to 
accurately describe the variations of the acoustic quantities in that 
direction. 

In Figs. 1-3, the distributions of calculational points 'n the surface 
of the body are presented. The arrows represent outward normals from the 
body as these plots represent a visual check of the geometric input data for 
the acoustic calculation programs. It will be noted here that if the driver 
plane, which doubles here as the duct exit plane as there is no straight duct 
section connected to the inlet lip, is non-dimensionalized as 1.0 then the 
nominal point spacings on tlie driver plane and the inner bellmouth are 0.04, 
0.025, and 0.02 respectively. The bellmouth geometry is terminated in all 
cases by a 2:1 ellipse as the theory requires a closed body. 
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The Computer Runs 


In Table II, a matrix of all of the computer r^ns done for the Langley 
Bellmouth is presented along with the five special runs done at the request 
or the contractor. This table also relates the cut-off ratios for the various 
modes run to the corresponding non-dimensional wave numbers at which the 
runs were performed. A standard plot of the sound pressure level calculated 
in the field is generated with each computer run. These plots are presented 
as Figs. 4-38. These figures are grouped according to the input mode used at 
the driver and are presented in ascending wave number order. Thus Figs. 4- 
12 are for an input mode of M-{. Figs. 13-20 for M 2, Figs. 21 26 for M .4, 
and Figs. 27-32 for M .8. All of these are done for the fifst radial mode 
input at the driver piane. The final set of these plots. Figs. 33-38 is for a 
tangential tnode of M 4 and the third radial mode input at the driver plane. 
For these plots, the data is calculated at 46 points in the field on a quarter 
circ'le having a radius of 20.36a. These points are evenly spaced from 0^, 
the centerline of the duct, to 90^ in 2° increments. Th.s is actually a rather 
course mesh when looking at higher radial modes as in Figs. 33-38 run for an 
input mode of (M,N)-(4,3) and this is why the valleys between the three 
lobed radiation pattern in the field are somewhat indistinct. Also, the points 
are calculated 20.36 radii from the duct entrance as this is where they are 
measured in Refs. 2 and 3 

If we now look at the "concept of angular location ot the 
predotninant acoustic radiation" v/e find that in genera!, this radiation peak 
moves towards the ("enterline of the duct for increasing wave numbers and 



away from tfie duct centerline for increasing tangential mode numbcis. 
These results, for all of tlie runs done for this inlet, are presented in Table 
III. For the modal input (4,3) the lobe that carries the maximum SPL in (dB) 
changes from tht; first (cocnted from the centerline of the duct) to the third 
at a cut-off ratio of 1.305. Thus, the bracketed values in Table III follow the 
first lobe even though it is no longer the predominant one (See Figs. 33-38.). 
It will be noted here again that since the acoustic quantities were only 
calcuated at 2^ increments from the centerline of the duct to 90^ that the 
angular values in Table III are probably only good to ^ 1 
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Reflection Coefficients 


Since there is no straight duct section connected to the inlet, the 
driver plane is the same as the duct exit plane. Therefore, the reflectii>ri 
coefficient defined as 

.1.1 f 1 ) 

I k ; .jn '* 

i •< 4’ + H. 

an 

where k is the wave number, is the acoustic [x->tential, 

and a I* IS the normal acoustic velocity (defined with an outw it d normal 

In 

to .he body), can be iwlculated directly from the results of the surface 
acoustic program without using the field program. Since the refle>-tion 
coefficient is complex, it is plotted as a modulus and a phase (in degrees) m 
two separate graphs with the same figure number. 

In Figs. 39 and 'rO, the reflection coefficient is plotted for an input 
mode of (1,1) for dll of the wave numbers (cut-off ratios) run. Some of the 
individual cases run are repeated in notti figures for ease of comparison. 
Similarly, in Figs. 41 and 42, the reflection coefficient at the duc t exit plane 
is plotted for an input mode of (2,1) at the driver plane for all of the cases 
run. Akso, Figs. 43 and 44 are for input modes of (4,1) and (8,1) respectively. 
For on input mode of (4,3) the results are plotted in six graphs designated as 
Fig. 43. Since the input mode has two .'ero crossings for the ai'ousti ■ 
poteiiti.il specified, the reflection coefficient has a discontinuity it these 
point.s. This is shown in Fig. 43 by artifically extending the data to the plot 
boundaries in the direction it was headed at these points. Since the trends 
were not .ilways in the saine dire.'tion at these points, the data was plotteil 
in separate graphs for clarity. 
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'n the next set of figures, the reflection coefficient is plotted for 
input modes of (1,1), (2,1), (4,1) and (8,1) for constant cut-off ratio. As 
before, each figure consists of two one for the modulus of the 

reflection coefficient and one for the phase. Figures 46-51 are for cut-off 
ratios of 1.005, 1.015, 1.064, 1.155, and 1.305 and 1.556 respectively. 

It wil* be noted here that the reason for the seemingly sirange 
behavior of the reflection coefficient in Figs. 44 and 51 for an input mode of 

(8.1) at a wave number of 15.01 (i.e. a cut-off ratio of 1.556) is that this is 
above the second radial mode cut-off wave number which is ka 14.12 for 

(8.2) . This is significant as the integral technique empluyed in the 

calculation of these values automatically considers all radial modes for the 
tangential mode specified. Thus, even though the pressure (i.e. potential) 
specified at the driver plane is for the (8,1) mode, the solution is picking up 
the effect of the (8,<^'' moo'e as well as for this wave number it is cut-on and 
has a cut-off ratio of 1.063, It is not possible to separate out this effect 
with this formulation of the problem as this would require the over 
specification of the boundary conditions on the driver plane. That is, it 
would require the specification of both the acoustic potential (i.e. pressure) 
and the normal acoustic velocity on the driver which cannot be done for an 
elliptic problem. 

This is also why the principle radiation peak shifts from the first to 
the third for an input mode of (4,3) (See Figs. 33-38 and Table III.). In this 
case, the (4,4) mode is cut-on at a wave number of 1 5.96 so that the last two 
wave numbers, 16.55 and 19.73, have this mode cut-on with cut-off ratios of 
1.037 and 1.236 respectively. 
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Radiation Patterns 


The radiation patterns are presented in Figs. ^=38 for the 3^ runs 
done for the .angely Bellmouth. Since in this form the results are not easy 
to compare, the relative SPL (Sound Pressure Level) in decibels and phase in 
degrees have been cross plotti.d so that trends can be seen rnore easily. The 
plots are of the relative values referenced in all cases to the point in die 
field where the maximum SPL occurs for each case. 

In Figs. 52 and 53, the results are plotted for all of the runs done for 
an input mode of M-(l,l). Each figure consists of three plots: tvvo for the 
relative SPL in dB and one for the phase. Note that the phase plot in Fig. 52 
has an expanded scale as compared to the phase plot in Fig. 53. This was 
done simply ^or the sake of clarity. This is also done in Fig. 54. As with 
the reflection coefficient plots some of the runs are repeated in Figs. 52 and 
53 for ease of comparison. Similarly, the results of all of the runs done for a 
modal input of M-(2,0 are presented in Figs. 54 and 55. Also, the results 
for input modes of (4,1) and (8,1) are presented in Figs. 56 and 57 
respectively. As mentioned before, the seemingly aberrant behavior of the 
results for a modal input of (8,1) at a non dimensional wave number of 
ka 15.01, clearly visible in Fig. 57b, is due to the fact that the (8,2) mode is 
cut-on at this frequency. This is one of the reasons why the results for an 
input mode of (4,3) for the last two wave numbers run, 16.55 and 19.73, are 
presented in separate plots in Fig. 58. Also, the valleys between the lobes in 
these plots would be more distinct if more points had been used in the field. 
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PirwUy, in Pigs. 59-6#, th« rwlativt results in the field are presented 
for modal Inputs of M>(1,1), (2,0, (#,1) and (S,l) Awhile the cut-off ratio is 
held constant as was done for the reflection coefficient plots. Again, t)>e 
relative SPL in dB is plotted In two ways In each figure. 
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Fig. 19 (M,N) = (2,1) 
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Fig. 25 (M,N) = (4,1) 
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Fig. 26 (M,N) = (4,1) 
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Fig. 27 (M,N) = (8,1) 
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Fig. 31 (M,N) = (8,1) 
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Fig. 36 (M,N) = (4,3) 
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Fig. 38 (M,N) = (4,3) 
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Fig. 62c 
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Fig. 63a 
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